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Abstract 

We introduce the Green’s functions technique as an alternative theory to the quantum regression theorem formalism 
for calculating the two-time correlation functions in open quantum systems. In particular, we investigate the potential 
of this theoretical approach by its application to compute the emission spectrum of a dissipative system composed 
by a single quantum dot inside of a semiconductor cavity. We also describe a simple algorithm based on the Green’s 
functions technique for calculating the emission spectrum of the quantum dot as well as of the cavity which can easily 
be implemented in any numerical linear algebra package. We find that the Green’s functions technique demonstrates 
a better accuracy and efficiency in the calculation of the emission spectrum and it allows to overcome the inherent 
theoretical difficulties associated to the direct application of the quantum regression theorem approach. 

Keywords: Quantum dot, semiconductor cavity, emission spectrum, Markovian master equation. Green’s functions 
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1. Introduction 

The measurement and control of light produced by 
quantum systems have been the focus of interest of 
cavity quantum electrodynamics [1, 2]. Specially, the 
emission of light powered by solid-state devices cou¬ 
pled to nanocavities is an extensive area of research due 
to its promising technological applications, such as in¬ 
frared and low-threshold lasers [3, 4], single and en¬ 
tangled photon sources [5, 6], as well as various appli¬ 
cations in quantum cryptography [7], and quantum in¬ 
formation [8]. Experiments with semiconductor quan¬ 
tum dots (QDs) embedded in microcavities have re¬ 
vealed a plethora of quantum effects and offer desirable 
properties for harnessing coherent quantum phenom¬ 
ena at the single photon level. Eor example, the Pur¬ 
cell enhancement [9], photon anti-bunching [10], vac¬ 
uum Rabi splitting [11] and strong light matter coupling 
[12]. These and many others quantum phenomena are 
being confirmed experimentally by observing the power 
spectral density of the light (PSD) emitted by quantum- 
dot-cavity systems (QD-Cavity). Thus, the PSD or so- 
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called emission spectrum is the only relevant informa¬ 
tion about the system which allows to study the prop¬ 
erties of light via measurements on correlations func¬ 
tions as stated by the Wiener-Khintchine theorem [13]. 
In order to compute the absorption or emission spec¬ 
trum in open quantum systems, more precisely, in QD- 
Cavity systems different approaches have been devel¬ 
oped from theoretical point of view. Eor example, the 
method of the thermodynamic Green functions which 
is applied to the determination of the susceptibilities 
and absorption spectra of atomic systems embbeded in 
nanocavities [14], and the time-resolved photolumines¬ 
cence approach whose application allows to determine 
the emission spectrum by consideration of an additional 
subsystem called the photon reservoir [15]. However, 
these methods have their own approximations and re¬ 
strictions and therefore are not widely used. Erequently, 
the emission spectrum in QD-Cavity systems is com¬ 
puted through the Quantum Regression Theorem (QRT) 
[16, 17, 18], since it relates the evolution of mean values 
of observables and the two-time correlation functions. 
It is worth mentioning that this approach can be diffi¬ 
cult to implement in a computer program, it due to that 
computational complexity of QRT approach increases 
significantly as the number of QDs or modes inside the 
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cavity, and the dimensionalities of the Hilbert spaces 
are large. In general, this approach is time-consuming 
due to that it requires to solve a large system of cou¬ 
pled differential equations, and numerical instabilities 
can arise. Moreover, theoretical complications can ap¬ 
pears related to dynamics of the operators involved, as 
we will point out in the next section. In spite of this, the 
QRT approach is widely used for theoretical works, for 
example, in studies of the luminescence spectra of cou¬ 
pled light-matter systems in microcavities in the pres¬ 
ence of a continuous and incoherent pumping [19, 20], 
and the relation between dynamical regimes and entan¬ 
glement in QD-Cavity systems [21, 22]. In the past, the 
Green’s functions technique (GFT) was successfully ap¬ 
plied for calculation of the micromaser spectrum [23], 
as a methodology in which the two-time correlation 
function is treated as a Green’s function that decays as 
the off-diagonal elements of the reduced density matrix 
of the system for a very specific initial conditions. Nev¬ 
ertheless, this approach has not been widely noticed in 
many significant situations in open quantum systems. 
Possibly, it is due to their work having a limitation of 
implementation. The purpose of this work is to present 
a simple, but efficient numerical method based on QRT 
formalism which overcome the inherent difficulties as¬ 
sociated to the direct application of the QRT, by solving 
the dynamics of the system in the frequency domain di¬ 
rectly. This paper is organized as follows: Section 2 re¬ 
view the theoretical background of quantum regression 
theorem and its relationship with the Green’s functions 
technique. Section 3 deal with a concrete application of 
our proposed method for computing the emission spec¬ 
tra of a dissipative QD-Cavity system. In section 4 we 
show the numerical calculations of the emission of spec¬ 
trum for the cavity and the quantum dot from both GFT 
and QRT approaches. Finally, we conclude in the last 
section. 

2. Theoretical background 

One of the most important measurements when the 
light excites resonantly a QD-Cavity system is the emis¬ 
sion spectrum of the system. From theoretical point of 
view, it is assumed that corresponds to a stationary and 
ergodic process which can be calculated as a PSD of 
light using the well-known Wiener-Khintchine theorem 
[13]. It states that the emission spectrum is given by the 
Fourier Transform of the correlation function (two-time 
expectation value) of the operator field c, 

1 r“ 

5(w)= -melim dr, (1) 
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where K{t) - {c\t + T)c(f)) and the normalizing fac¬ 
tor is the population tic at the steady-state. In or¬ 
der to calculate the two-time expectation value is fre¬ 
quently used the QRT which states that if a set of 
operators {Oj{t + t )] satisfy the dynamical equations 
i{Oi{t+r)) = 2,. L,/GJf+T)> then ^<G,(f+T)G(f)) = 
Yjj Lij{Oj(t + T)O(t)) is valid for any operator (9(f) at ar¬ 
bitrary time f. It is worth mentioning that vality of this 
theorem holds whenever a closed set of operators are in¬ 
volved in the dynamics. In general, to obtain the closed 
set of operators can be difficult or an impossible task, 
since it must be added as many operators as necessary 
in order to close the dynamics of the system. For ex¬ 
ample, in order to compute the emission of spectrum in 
a simple model of QD-Cavity system [20, 21] two new 
operadors are required due to that the field operators in 
the interaction picture does not lead to a complete set. 
Before we consider the Green’s functions technique, we 
will briefly describe the calculation of the QRT in an 
alternate form which will be the starting point in the 
following section. Lets consider a system operador A 
which does not operate on the reservoir, then its single¬ 
time expectation value in the Heisenberg picture is given 

by 

<A(f + t)) = T rs®«[A(f H- T)pj««(f)]. (2) 

The operator ps»R(t) - ps (f) ® puit) depics the com¬ 
posite density operador of the system and reservoir. 
It is worth pointing out that the dynamics of the sys¬ 
tem depends directly on pj®«(f) for all times, but the 
validity of the Markovian approximation requires that 
the state of the system is sufficiently well described by 
Ps (f) = TrR(pse>R(t)), therefore it is sufficient to write 
Ps®«(0 = Ps(t) ® P«(f)- In what follows, we change 
to the Schrodinger representation using A(f H- r) = 
U’'(t+T, t)A(t)U{t + T, t) with U(t + T, t) being the unitary 
time-evolution operator, and after tracing over degrees 
of freedom of the reservoirs, we have 

{A{t + T))^Trsm)ps{t + T)], ( 3 ) 

where the reduced density operador for the system is 
given by ps(f -H t) = Tr«[0(f -H t, f)pj®^(f)0'(f H- r, f)]. 
Then, if the pj(f + t) satisfies the Markovian master 
equation dpsit + t) I dr = X.ps(t + r) with X. the Li- 
ouvillian superoperator, the evolution of (A(f -H t)) can 
be computed by solving the dynamics of the master 
equation. To calculate the two-time correlation func¬ 
tion (A(f -I- r)B(f)) where A(f -i- r) and B{t) are arbitrary 
Heisenberg operators, we proceed in a similar manner, 
it is, 

<A(f -H T)B(f)) = Trsm[Mt + T)B(f)Ps®R(0], 

= Trs[Ait)Git + T)], (4) 
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where we have used the well-known composition and 
inversion properties of the evolution operator. Then, the 
two-time operator is given by 

G(t + t)^T rR[U(t + T, t)B{t)psm{t)UHt + t, t)]. (5) 

By comparison of the Eq. (3) and Eq. (4), we find that 
G(t + t) is an operator that obeys the same dynamical 
equations as ps (t + r), but as function of t. It is, dG(t + 
T)/dT - £.G{t + t) with the boundary condition G(f) = 
B(t)ps(t) at arbitrary time t. Hence, in the long-time 
limit the QRT reads, 

lim<A(f + T).B(f)) = Trs [AG(t)] (6) 

f—»oo 

where G(t) = 7’rR[G(T)Bp^^^^G^(T)] is the Green’s 
functions operator, and the operators A, B and p\Jr are 
written in the Schrodinger representation. The super¬ 
script ”(ss)” refers to the steady state of the reduced 
density operator of the system. After taking the Laplace 
transform on Eq. (6), we obtain an expression for the 
emission of spectrum in terms of the Green’s functions 
operator, it is, 

S{co) = — me Trs [AG(/w)]. (7) 

mic 

Prior to leaving this section, we mention that this equa¬ 
tion will be considered for computing the emission 
spectrum due to the cavity as well as the quantum dot, 
e.g. by considering the photon and fermionic operators 
in a separated way. Therefore, we will describe in the 
next subsetion a general approach that can be applied 
for both cases. 

2.7. Algorithm for the Green’s functions technique 
Before to describe a simple algorithm for calculating 
the emission spectrum, we take into account that the dy¬ 
namics for both opertors G(t) and ps (t) are governed by 
the same Master equation, i.e., dG(T)ldT - J1 G{t) with 
£. the Liouvillian superoperator, that efectivelly has a 
larger tensor rank than the reduced density operator of 
the system. So, we can write the dynamical equations 
for the Green’s functions operator in a component form: 

( 8 ) 

together with the initial condition Gjj{Q)- The symbol 
(f is a composite index for labeling the states of the re¬ 
duced density operator of the system, e.g. for indexing 
both matter and photon states in the QD-Cavity system, 
see section 3 for example. Hence, G^ and acts as 


a column vector and a matrix in this notation, (ii) To 
obtain the solution to the Eq. (8) in frequency domain 
via the Laplace Transform, it is -Gq.(0) = — 

ico5^p)Gp{ico), (iii) We perform the invertion of the ma¬ 
trix - (iatS^p - X.^p) and finally, the spectrum of 
emission is computed in terms of the initial conditions 
given by, 

Gp(iaj)^^M-/,GdO). (9) 

a 

These initial conditions are easily obtained by evaluat¬ 
ing the Green’s function operator at t = 0. 

3. Application to the quantum dot-cavity system 

3.1. Model 

In order to apply our proposed method for calculat¬ 
ing the emission spectrum in QD-Cavity system, we 
will consider a simple but illustrative system composed 
of a quantum dot interacting with a confined mode of 
the electromagnetic field inside a semiconductor cavity. 
This quantum system is well described by the Jaynes- 
Cummings Hamiltonian [24] 

Hs — a)x(T^O' + {(^x - A)a^d + gip-a^ + a&^), ( 10 ) 

where the quantum dot is described as a fermionic sys¬ 
tem with only two possible states, |G) and |A) are the 
ground and excited state, respectively. & - |G)(A| and 
d ((T* = |2f)(G| and aQ are the annihilation (creation) 
operators for the fermionic system and the cavity mode. 
g is the light-matter coupling constant, and we have set 
h-\. We also define the detuning between frequencies 
of the quantum dot and the cavity mode as A = cjx-cja, 
moreover a>x is the energy to create an exciton and 0 )^ 
is the energy associated to the photons inside de cavity, 
respectively. This Hamiltonian system is far away for 
describing any real physical situation since it is com¬ 
pletely integrable [25] and no measurements could be 
done since the light remains always inside the cavity. 

In order to include the effects of environment on the dy¬ 
namics of the system, we consider the usual approach to 
model an open quantum system by considering a whole 
system-reservoir hamiltonian which is frequently split- 
ted in three parts. The first part corresponds to the sys¬ 
tem of quantum dot-microcavity. The second part is the 
hamiltonian of the reservoirs and finally, the third part 
which is a bilinear coupling between the system and the 
reservoirs [26]. After tracing out the degrees of free¬ 
dom of all the reservoirs and assuming the validity of 
the Born-Markov approximation, one arrives to a master 
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equation for the reduced density matrix of the system. 


dps 

dt 


- i [ps ,Hs] + -(2apsa^ - a^aps - psd^d) 

y f f t 

+ -{2&PS&' -a-'&ps -psd-'d-) 

p 

+ - d-a^ps - Psd-&'^)- ( 11 ) 


Where y is the decay rate due to the spontaneous emis¬ 
sion, K is the decay rate of the cavity photons across the 
cavity mirrors, and P is the rate at which the excitons 
are being pumped. Fig. 1 shows a scheme of the sim¬ 
plified model of the QD-cavity system showing the pro¬ 
cesses of continuous pumping P and cavity loses k. The 
physical process begin when the light from the pump¬ 
ing laser enters into the cavity and excites one of the 
quantum dots in the QD layer. Thus, light from this 
source couples to the cavity and a fraction of photons 
escapes through the partly transparent mirror from the 
cavity and goes to the spectrometer for measurements 
of the emission of spectrum. 

A general approach for solving the dynamics of the cou¬ 
pled system, consist of writting the Bloch equations for 
the reduced density matrix of the system in the bared 
basis. It is, an extended Hilbert space formed by taking 
the tensor product of the state vectors for each of the 
system components, {|G), |X)) ® {|n))J^Q. In this basis, 
the reduced density matrix ps can be written in terms of 
its matrix elements as psanfim - {an\ps\Pni). Hence, the 
Eq. (11) explicitly reads. 


dpSan,pm 

dr 

i[(aJx - A)(m - n)psan,pm 

+ 

tOx{dpxPSan,Xm - daXPSXnfim)\ 

+ 

Vm H- \5pxPSan,Cm+l 

+ 

'^dpcPSan.Xm-l) 

- 

( ^daCPSXn-lfim 

+ 

V« -H l5axPSCn+\pn^ 

+ 

-^2 yjim + 1)(« + l)psan+\fim+\ 

- 

ill + fyi)pSan.fiin^ ~ f^idaXPSXn,^}}! 

- 

'^dacdpcPSXn.Xm + dpxPSan,Xm) 

+ 

^{^daxdpxPSCn.Cm - daCPSGnfim 

- 

dpGPSan,Grri^‘ (1^) 


Note that we use the convention that all indices writ¬ 
ten in greek alphabet are used for matter states and take 


values |G), |X), and the indices written in latin alpha¬ 
bet are used for Fock states and take values 0,1,2,3 .... 
Additionally, it is worth mentioning that our proposed 
method does not require to solve a system of coupled 
differential equations, instead of it, we solve a reduced 
set of algebraic equations that speed up the numerical 
solution. 

Prior to leaving this section, we point out that the num¬ 
ber of excitations of the system is defined by the opera¬ 
tor A = a^a+a-^&. The closed system and the number of 
excitations of the system is conserved, i.e., A] = 0. 
It allows us to organize the states of the system through 
the number of excitations criterion such that the density 

matrix elements pGn,G«.Pxn-i,Xn-i Pcn.Xn-i andpx«-i,G« 

are related by having the same number of quanta. It 
is, subspaces of a fixed number of excitation evolve 
independently from each other. The Fig. 2 shows a 
schematic representation of the action of the dissipa¬ 
tive processes involved in the dynamics of the system 
according to the excitation number (Nexc)- 


Substrate 



Figure 1: Scheme of the simplified model of the QD-cavity system 
showing the processes of continuous pumping P and cavity loses k. 


3.2. Emission spectrum of the cavity 
In order to compute the emission spectrum of the cav¬ 
ity, we will consider the two-time correlation function 
according to the Eq. (6) for the photon operator as fol¬ 
lows: 

K(t) - lim(a^(f + T)a{t)). (13) 

t—^oo 

After performing the partial trace over the degrees of 
freedom of the system we have that 

K(t) = sj{l + l)(m + l)Trfi[UaijsmiT) 

X + (14) 
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Figure 2: Ladder of bared states for a two level quantum dot coupled 
to a single cavity mode. The double headed green arrow depics the 
matter coupling constant g, dashed red lines the emission of the cavity 
mode K, solid black lines the exciton pumping rate P and solid blue 
lines the spontaneous emission rate y. 


where the matrix elements for the time evolution op¬ 
erator are given by Uaifimij) - {al\U{T)\pm) and 
= {'Yn\U"^{T)\al + 1). In what follows, we 
assume the validity of the Markovian approximation, it 
means that the correlations between the system and the 
reservoir must be unimportant even at the steady state. 
Thus, the density operator system-reservoir can written 
P^s»R ~ ® Pr^^ which implies that 


In order to identify these matrix elements, it should con¬ 
sider that for the QD-Cavity system the dynamics of 
the all coherences asymptotically vanish and remains 
only the reduced density matrix elements which are 
ruled by the number of excitations criterion, i.e. pcn,cn, 
Pxn-\,Xn-\, PGn,Xn-\, Pxn-i,Cn- Then, the Eq. (15) can be 
written as follows, 

</3m H- llp^^^^lyn) = Ps^^(dpc5yG^m+l,n 

+ 5px5yx5m,n-\ + SfjcSyX^m^n 
+ ^0xSycSm+l,n-\)p%l+i^y„- 
( 18 ) 

By replacing the Eq. (18) into Eq. (17) we hnd that the 
Green’s functions operator explicitly reads 

G(t) = TrR[u(T)p^^'^ Yj + ^{\Gm}{Gm + 1| 

m 

® PsCm+l,Cm+l |Xm)(Xm -H ^P^^Xm+\ ,Xm+l 

-H \Gm){Xm\p^^2m+\,Xm + |Xm><Gm H- 2| 

>< PfL.,p,n.2)uHT)\ (19) 

Note that from this expression is easy to identify the 
nonzero matrix elements of the Green’s functions oper¬ 
ator that contribute to the emission spectrum. Einally, 
after performing the Laplace transform we have that the 
emission spectrum of the cavity is given by 


(fim + llp^-'^’^lyn) = p^^"^{j3m + Ilp^-^'V”)- (15) 

Replacing the previous expression in Eq. (14), it is 
straightforward to shows that the two-time correlation 
function reads 

K(t) = Yj ViTT<al|G(T)|al H- 1), (16) 

al 

where the Green’s functions operator G(t) is given by 
G(t) = 7’rfi[G(T)p^^^^ ^ (VmTTlySmXynl 

/3ymn 

X {fim+l\pf'>\yn))u\T)]. (17) 

As we pointed out in section 2, this operator must obey 
the same master equation as the reduced density oper¬ 
ator of the system. In fact, the terms that only con¬ 
tribute in the Eq. (16) are given by the matrix elements 
Gpm,yn{T) = (j3mlG(T)lyn} of the Green’s functions oper¬ 
ator. This is due to the fact that the projection operator 
IjSmXynl enter into G(t) in the same way as into the re¬ 
duced density operator of the system. 


S(a>) - -^ V/ -H 1 (Gc/p/+i(ioj) 

nur ^ ^ 

I 

+ Gxi,xi+i{i(^) + GcuxiiitG) 

+ Gxi,ci+2(i<^))- (20) 

It is worth mentioning that the initial conditions may be 
obtained by evaluating the Green’s function operator at 
T = 0, then using the fact that the time evolution opera¬ 
tors become the identity and T rR[p^^^^] = 1, we obtain a 
set of initial conditions given by 


Ggi,ci+i(0) - 

Vi -r 


Gxi,xi+i(0) = 

Vi -r ^f' 5 x/+l,X/+l’ 


GGijiiiO) = 

Vi -r 


Gx/,G/+2(0) = 

Vi -r if' 5 x/+l,G/+ 2 ‘ 

(21) 


Note that this set of initial conditions corresponds to the 
asymptotic solution of the Bloch equations for the re¬ 
duced density matrix of the system. 
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3.3. Emission spectrum of the quantum dot 

In order to compute the emission spectrum of the quan¬ 
tum dot, we will consider the two-time correlation func¬ 
tion given by Eq. (6), but for the case of the matter op¬ 
erator; 

K{t) - lim((3-^(f + T)&(t)). 

t~-*oc 

It is straightforward to show after performing the partial 
trace over the degrees of freedom of the system that the 
two-time correlation function reads 

KiT)^Y,^,x{Gl\GiT)\al}, (22) 

al 

where the Green’s functions operator G(t) is given by 
G(t) = TrR[u(T)Y^(^dfix\Gm}{ym\ 

^ymn 

X (J3m\pfJ^\yn))uHT)]. (23) 

Assuming again the validity of the Markovian approx¬ 
imation and taking into account the number of excita¬ 
tions criterion, we have that the density operator system- 
reservoir can be written as: 

+ + ^l3C^yxSm,n+\ 

+ dpxdycdm,n-'i)p^sp^yn- 

(24) 

By inserting the Eq. (24) into Eq. (23) we find that the 
Green’s functions operator explicitly reads 

G(r) = 

m 

+ |Gm)<Gm+l|p(->„_^^^,)Gt(T)]. (25) 

Analogously as in section 3.2, we identify the nonzero 
matrix elements of the Green’s functions operator that 
contribute to the emission spectrum and after perform¬ 
ing the Laplace transform the emission spectrum of the 
quantum dot is given by 

S (w) = -— ^ (Cci.xiiit^) + Gci.ci+iiitf))- (26) 

where Ug- = is the normalizing factor at the 

steady-state. Taking into account that the initial con¬ 
ditions are obtained by evaluating the Green’s function 
operator at t = 0, we have the time evolution operators 
become the identity and = 1, thus, we obtain 


a set of initial conditions given by 
GcuxiiO) = P^sxuxp 

Gg/,g/+i(0) = Psxi.ci+v 

Gxz,z/+i(0) = 0, 

Ggz+2.x;(0) = 0. (27) 

4. Results and Discussion 

In this section, we compare the numerical calcula¬ 
tions based on GET and QRT approach for the emis¬ 
sion spectrum of the cavity as well as the quantum dot. 
Due to that the QD-Cavity system can display two dif¬ 
ferent dynamical regimes by changing the values of the 
free parameters of the system and transitions between 
these two regimes can be achieved when the loss and 
pump rates are modihed. Particularly, in the strong cou¬ 
pling regime the relation Pjx ^ g holds and the re¬ 
lation PjK » g remains valid for the weak coupling 
regime. Fig. 3 shows the numerical calculations of the 
emission spectrum due to the cavity in the weak cou¬ 
pling regime, the parameters values are g = 1 meV, 
y — 0.005 meV, k — 0.2 meV, P = 0.3 meV, A — 2 meV, 
cJa - 1000 meV. Panel (a) shows the emission spec¬ 
trum for the GET compared to the QRT approach. Panel 
(b) shows the quantity \S(a>)GFT - S{(jj)qrt\ as a mea¬ 
surement of the error between the numerical calcula¬ 
tions in the emission spectrum. For this set of param¬ 
eters values, we can easily to identify two peaks asso¬ 
ciated to the modes of the cavity and the quantum dot, 
it is Ota ~ 998.3mey and a>x ~ lQ00.3meV, respec¬ 
tively. Fig. 4 shows the same calculations as in Fig. 3, 
but in the strong coupling regime and the parameters 
values are g = ImeV, y = 0.005 meV, k — 2meV, 
P - 0.005 meV, A = O.OmeV, oJa - lOOOmeV. In the 
case of resonance, the modes associated to the cavity 
and the quantum dot do not match, but repel each other, 
resulting in a structure of two separate peaks a distance 
2g X 2meV. Fig. 5 shows the numerical calculations of 
the emission of spectrum due to the quantum dot with a 
high value of the rate k - 5 meV and a smaller, although 
non negligible, pumping P = 1 meV. The rest of param¬ 
eters values are g = 1 meV, y - 0.1 meV, A = 5 meV, 
CJa - 1000 me V. 

We observed that our numerical method based on GET 
is in full agreement with the QRT approach and repro¬ 
duces very well the spectrum of emission associated 
with this system. The quantity |5 (cu)gf7- - S(cj)qrt\ 
shows the discrepancy between both methods which is 
the order of 10~^ - 10“^ as it is seen in Fig. 4 and Fig. 3. 
Mainly, the discrepancy between both methods is due to 
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the numerical errors accumulated in the numerical inte¬ 
gration of the Bloch equations in the QRT approach, it 
causes some differences in the spectrum with respect to 
the results computed by the GFT. Note that, there is no 
integration of any equations in the GFT, therefore, we 
expect a more accurate spectrum of emission. We men¬ 
tion that for the numerical calculations based on QRT 
approach, we have followed the Ref. [21]. In order to 
test the performance of the numerical method, we re¬ 
gard four calculation times for computing the emission 
spectrum of the cavity based on GFT and QRT approach 
at different excitation numbers. Table 1 shows in hrst 
column the excitation number, i.e. the truncation level 
in the bare-state basis for the numerical calculations in¬ 
volved. Second and third column show the results of 
elapsed time in seconds for both the GFT and QRT ap¬ 
proach, respectively. Note that for comparison purposes 
all numerical calculations were performed at the same 
truncation level i.e. Nexc - 10. Additionally, we have 
solved numerically the Bloch equations (see Eq. (12)) 
until time tmax - 2'^ in order to obtain a good resolu¬ 
tion in the frequency domain for the QRT approach, i.e. 
Aw =» 0.048. Hence, we have evaluated the emission 
spectrum for the GFT in a grid with the same resolution 
in the frequency domain (we emphasize that QRT ap¬ 
proach is time consuming due to the number of coupled 
differential equations to be solved, rather than the num¬ 
ber of evaluations in the grid size used). In addition, the 
numerical calculations were catTied out with the same 
parameters values as in Fig. 4 for both GFT and QRT 
approach. We found that our numerical approach based 
on GFT is very efficient and accurate for calculating the 
emission spectrum in QD-Cavity systems. Moreover, 
this method can easily be implemented in the numerical 
linear algebra packages as well as in any programming 
language. 


Table 1: Comparison of calculation times between the Green’s Func¬ 
tions Technique (GFT) and the Quantum Regression Theorem (QRT) 
in the numerical calculation of the emission spectrum of the cavity. 
The calculations were made using a commercial Intel(R) Core(TM) 
il - 4770 processor of 3.4 GHz x8, and 12 GB RAM. 


Excitation 

number 

Calculation times 
for the GET (s) 

Calculation times 
for the QRT (s) 

5 

0.4 

92.5 

10 

2.0 

273.4 

20 

14.0 

390.2 

40 

100.2 

673.8 



Figure 3: Panel (a) shows a comparison of the emission spectrum of 
the cavity. The numerical calculation based on the Green’s functions 
technique (GFT) is shown as solid blue line and the con'esponding 
numerical calculation based on quantum regression theorem (QRT) 
approach is shown as dashed red line. Panel (b) shows in solid blue 
line the quantity \S{(jj)gft ~ S(cS)qrt\ as a measure of the dilference 
in the numerical calculations of the emission spectrum of the cavity 
between two methods. 


5. Conclusions 

We have developed the Green’s function technique as 
an alternative methodology to the QRT for calculating 
the two-time correlation functions in open quantum sys¬ 
tems. In particular, we have shown the performance of 
the Green’s function technique by calculating the emis¬ 
sion spectrum in an open quantum system composed by 
a quantum dot embedded in a microcavity. This theoret¬ 
ical approach is rather general and allows to overcome 
the inherent theoretical difficulties presented in the di¬ 
rect application of the QRT, i.e., to hnd a closure con¬ 
dition on the set of operators involved in the dynamics 
equations, by considering that all coherences asymptot¬ 
ically vanish, and remains only the reduced density ma¬ 
trix elements which are ruled by the number of excita¬ 
tions criterion. We have shown that the Green’s func¬ 
tion technique offers several computational advantages, 
namely, the speeding up numerical computations via a 
transformation of the dynamics of the master equation 
in a set of linear algebraic equations, which are effi¬ 
ciently solvable by a numerical linear algebra routine, a 
faster convergence and signihcant reduction of compu- 
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Figure 4: Panel (a) shows a comparison of the emission spectrum of 
the cavity. The numerical calculation based on the Green’s functions 
technique (GFT) is shown as solid blue line and the coiTesponding 
numerical calculation based on quantum regression theorem (QRT) 
approach is shown as dashed red line. Panel (b) shows in solid blue 
line the quantity \S{co)QfT - S{cj)qrt\ as a measure of the difference 
in the numerical calculations of the emission spectrum of the cavity 
between two methods. 

tational time since the emission spectrum is calculated 
as a sum of terms of non-diagonal matrix elements of 
the reduced density operator of the system. We men¬ 
tion that our methodology can be extended for calcu¬ 
lating the emission spectrum in significant situations of 
quantum dots in biexcitonic regime or involving cou¬ 
pled photonic cavities. 

6. Acknowledgements 

This work was financed by Vicerrectoria de investiga- 
ciones of the Universidad del Quindio within the project 
with code 659, and by Colciencias within the project 
with code 110156933525, contract number 026 - 2013, 
and HERMES code 17432. 

References 

[1] H. Walther, et. al. Rep. Prog. Phys., 69 (2006) 395603. 

[2] A. Kavokin, et. al. Oxford University Press (2007) 

[3] H. Altug, D. England and J. Vuckovic Nature Phys. 
2 (2006) 484. 

[4] Y. Mu and C. M. Savage, Phys. Rev. A, 46 (1992) 5944 

[5] R. M. Stevenson, et. al. Nature 179 (2006) 439 


Figure 5: Panel (a) shows a comparison of the emission spectrum of 
the quantum dot. See text and caption of Fig. 4 for details. 


[6] T. M. Stace, G. J. Milburn, and C. H. W. Barnes, Phys. Rev. 
B, 67 (2003)085317 

[7] N. Gisin, et. al.. Rev. Mod. Phys. 74 (2002) 145 

[8] C. Monroe, Nature, 146 (2002) 238 

[9] Y. Todorov, et. al., Phys. Rev. Lett., 99 (2007) 223603 

[10] J. Wiersig, et .al Nature 460 (2009) 245. 

[11] G. Khitrova, et. al. Nat. Phys., 2 (2006) 81. 

[12] J. P. Reithmaier, et. al. Nature, 432 (2004) 197. 

[13] Mandel L. and Wolf E., Optical Coherence and Quantum Optics, 
Cambridge University Press, (1997) 

[14] O. Jedrkiewicz and R. Loudon, J. Opt. B: Quantum Semiclass. 
Opt., 2 (2000) R47. 

[15] N. V. Hieu and N. B. Ha, Adv. Nat. Sci: Nanosci. Nanotech- 
nol., 1 (2010)045001. 

[16] Walls D. F. and Milburn G. J. 1994 Quantum Optics (Berlin: 
Springer - Verlag). 

[17] M. Lax, Phys. Rev. 145 (1966) 110. 

[18] S Swain, J. Phys. A: Math. Gen., 14 (1981) 2577. 

[19] E. del Valle, F. P. Laussy, and C. Tejedor, Phys. Rev 
B, 79 (2009) 235326. 

[20] N. Quesada, H. Vinck-Posada and B. A. Rodriguez, J. 
Phys.: Condens. Matter, 23 (2011) 025301. 

[21] C. A. Vera, N. Quesada, H. Vinck-Posada and B. A. Rodriguez, 
J. Phys.: Condens. Matter, 21 (2009) 395603. 

[22] N. Ishida, et. al. Sci. Rep., 3 (2013) 1180. 

[23] T. Quang, et. al., Phys. Rev. A, 48 (1993) 803 

[24] E. T. Jaynes and F. W. Cummings, Proc. IEEE, 51 (1963) 89. 

[25] M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge: 
Cambridge University Press (1996) 

[26] J. I. Perea, D. Porras and C. Tejedor,Phys. Rev. 
B, 70 (2004) 115304. 


8 












